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Cementite Fe3C is studied with first-principles calculations. 
Two possible positions of carbon atoms in the iron sublat- 
tice are considered: with prismatic or octahedral environment. 
Mossbauer spectra (MS) with parameters calculated for both 
modifications are simulated above and below the Curie temper- 
ature. A possibility to detect the change in carbon position upon 
annealing from MS is discussed. It is shown that this is hardly 
possible using a standard approach to treatment of MS, but it 
can be seen in more subtle details of the MS below the Curie 
temperature, such as widths and positions of separate lines. 
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The properties of cementite, i.e. iron carbide Fe%C, are 
being studied for a long time. Originally, cementite was in- 
vestigated in order to improve the mechanical properties 
of steels. Its magnetic properties have received some atten- 
tion in connection with nondestructive methods of the steel 
control. A recent interest in the Fe^C properties is aroused 
by new possibilities of obtaining metastable compounds by 
mechanical alloying, implantation and so forth (see, for ex- 
ample, OEEIE)). A R additional impact has been given by 
the studies dealing with the chemical composition of the 
Earth's core |5]|6). Besides, the monophase Fe^C by itself 
has attracted considerable interest due to some peculiari- 
ties in the bulk-modulus behaviour Q and the instability 
of the magnetic state under pressure [09). 

The structure of cementite is believed to be known well 
from the experiments on the X-ray diffraction, for exam- 
ple, iflOl . But many magnetic, Mossbauer and EELS ex- 
perimental data [ i 1112:131, conducted for non-equilibrium 
samples, lead us to think that its structure varies depending 
on the mechanical and thermal treatment. As these varia- 



tions are not detected by the X-ray diffraction, one can sup- 
pose that there are cementite modifications differing in the 
carbon positions, which do not noticeably affect the X-ray 
diffraction patterns because of the large difference in the 
scattering factors between iron and carbon. The neutron 
diffraction experiments conducted recently on the well- 
annealed samples iFBl showed that carbon atoms occupy 
the prismatic positions. So, after annealing at high temper- 
atures (above 700 K IU0II14|[TD ) cementite is in the ground 
state and has a sole, well-determined, structure. The struc- 
ture of the non-equilibrium cementite comes into question. 

According to the structural considerations in [15], car- 
bon can occupy four positions between the iron sites, that 
are called prismatic and octahedral, and distorted prismatic 
and octahedral. The carbon in the two last positions is 
closer to iron atoms and these modifications appear to be 
improbable. Fig.[T]displays the structure of cementite, with 
two possible carbon positions shown together. 

In this work, first-principles calculations have been 
conducted for cementite with different positions of carbon 
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Figure 1 Double unit cell of cementite. Large black and 
grey balls show iron Fel (S position) and Fell (G position) 
atoms; small white or dark grey balls represent carbon in 
prismatic or octahedral environment, correspondingly. 



in the lattice. Earlier in JT6], using such calculations, we 
have received that the prismatic modification of cemen- 
tite has a smallest energy and is its ground state, and ex- 
plained an unusual behavior of the coercive force in the 
mechanically alloyed cementite on the annealing temper- 
ature as a transition of cementite from octahedral mod- 
ification with small magneto-crystalline energy Ema to 
the prismatic one with large -Ema- At the same time, we 
do not know unambiguous direct experiments demonstrat- 
ing existence of the octahedral modification of cementite. 
The X-ray diffraction experiments do not show the car- 
bon position with a necessary accuracy, and the neutron 
powder diffraction experiments have been done for well- 
annealed samples being in the ground state without traces 
of metastable phases. So, the aim of the paper is to calcu- 
late and compare the parameters of hyperfine interaction of 
cementite with different position of carbon atoms, to sim- 
ulate the Mossbauer spectra (MS) and to estimate the pos- 
sibility of experimental detection of the octahedral modifi- 
cation of cementite with the help of the Mossbauer spec- 
troscopy. 

Some theoretical investigations of cementite have been 
already carried out (see, for example, (6l[8l rT7lfT8ll ). Only 
in ifTTl an attempt to obtain and compare some character- 
istics of the prismatic and octahedral cementite modifica- 
tions has been done, but no analysis on experimental de- 
tection has been conducted. 

The crystal structure of Fe^C is an orthorhombic lat- 
tice, we use the Pbnm space group. The unit cell con- 
tains 12 iron atoms and 4 carbon atoms (a double unit cell, 
drawn with the help of the package |[T9l in Fig. [Tj shows 
carbon positions in both prismatic and octahedral environ- 
ment). 



The calculations presented in this paper have been con- 
ducted by the full-potential linearized augmented plane 
wave (FLAPW) method in the WIEN2k package 1201 . 
The generalized gradient approximation of the exchange- 
correlation potential (GGA) |2T] is used in this work. The 
details of the calculational procedure can be seen in lfl6l . 
Calculated are the Fe atomic magnetic moments, hyper- 
fine magnetic fields (HFF), isomer shifts (IS), the electric 
field gradient (EFG) at the Fe nuclei. V zz and V^ x are the 
maximum- and the minimum-in-magnitude components of 
the tensor of the EFG; i] = (V xx — Vy y )/V zz ; 6 and <fr are 
the polar and azimuthal angles of the HFF in the local coor- 
dinate system of the EFG. Using the parameters calculated, 
we have simulated the Mossbauer spectra for polycrystals 
(which is usually characteristic of cementite) with the help 
of the program described in 11221 . 

Table[TJgives data (cell parameters, magnetic moments, 
Mossbauer parameters) for each system. A difference be- 
tween the magnetic moments in the prismatic and octa- 
hedral cementite is found numerically small, so magnetic 
measurements are not expected to determine which carbon 
position is realized in samples. 

A difference in the EFG between the two modifica- 
tions draws attention. One could expect that this differ- 
ence should manifest itself in the MS above the Curie tem- 
perature. The simulation of the spectra with zero hyper- 
fine field and other parameters from Table 1 (see Fig. |2|i 
shows, however, that practically no difference can be seen 
between the two modifications. This occurs because the 

Table 1 The calculated structural parameters of the unit 
cells - lattice parameters a, b and c (nm) and positions 
(fractional coordinates) of the atoms, spin M sp i n and or- 
bital M oro magnetic moments (he), the Fermi-contact 
H sp i n and orbital H or b contributions to the hyperfine mag- 
netic field H t ot = H S pi„ + H oro (T), isomer shift IS 
(mm/s), electric field gradient EFG V zz (V/m 2 ), asym- 
metry parameter i], polar and azimuthal angle between 
H to t and main axis of the EFG tensor. 





prismatic 


octahedral 


a, b, c 


0.4490, 0.5047, 0.6743 


0.4711, 0.5111, 0.6960 


number 


Fe I, 


Fell, 


Fel, 


Fe II, 


of atoms 


4 atoms 


8 atoms 


4 atoms 


8 atoms 


Jospin 


1.97 


1.91 


2.09 


1.90 


M orb 


0.05 


0.04 


0.05 


0.03 


Hspin 


-25.1 


-24.7 


-24.3 


-23.7 


H or b 


3.1 


2.7 


3.1 


1.9 


Htot 


-22.0 


-22.0 


-21.2 


-21.8 


IS 


0.16 


0.17 


0.22 


0.14 


v„ 


3.07 


1.37 


1.99 


-2.09 


V 


0.05 


0.83 


0.996 


0.23 


<t> 


90 


138.6 


90 


73.1 


e 





11.2 


90 


26.6 
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Figure 2 The simulated 
Mossbauer spectra of ce- 
mentite with parameters 
from Table 1 above the 
Curie temperature (the 
HFF's equal zero). Curves 
1 are the subspectra of two 
positions of iron in the 
octahedral modification, 
curves 2 are for the prismatic 
modification, curves 3 are 
for the summary spectra of 
the octahedral or prismatic 
modifications. 
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Figure 3 The simulated 
Mossbauer spectra below 
the Curie temperature. 
Designations are as in Fig. 2 



difference wears off due to summation of two subspectra 
from nonequivalent positions Fe I and Fe II. 

Below the Curie temperature (see Fig. |3), the spec- 
tra can be distinguished better, as the difference becomes 
clearer in the outermost lines. The subspectra are essen- 
tially affected by the combined hyperfine magnetic and 



quadrupole interaction. The summary spectra cannot be 
also decomposed at a glance, but some traces of a com- 
bined interaction still exist and manifest themselves in the 
distances between the peaks of the simulated MS (see Ta- 
ble We remind that in the absence of the quadrupole 
interaction the distances between the first and second lines 
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and the fifth and sixth lines in the Mossbauer sextet are 
equal. 

Table 2 The distances between the first and second (R12) 
and the fifth and sixth (R^e) lines of the Mossbauer spec- 
trum in the case of the combined hyperfine interaction for 
the positions Fe I and Fe II and for the summary spectra. 





prismatic 


octahedral 




Fel 


Fe II 


summary 


Fel 


Fell 


summary 


R12 


1.71 


1.43 


1.52 


1.43 


1.33 


1.38 


Rse 


1.21 


1.52 


1.50 


1.43 


1.64 


1.57 



In an actual non-equilibrium sample, an even more 
complicated case should be realized, as there should be 
a mixture of the Fe^C modifications, and only a scrupu- 
lous mathematical processing of spectra with a preliminary 
model of the relations between the HF parameters obtained 
theoretically for two modifications can give a necessary ev- 
idence of the octahedral modification. 

The quadrupole splitting above the Curie temperature 
and its absence below have been also observed earlier in 
experiments of [23 1. The parameter of quadrupole splitting 
measured above Tc in 11231 corresponds to our calculation 
(see lfl6l ). However, the authors supposed that the disap- 
pearance of the quadrupole interaction in the spectrum be- 
low Tc is a result of a special angle between EFG and HFF 
(55 degree). Our calculations show that this is not true for 
all positions of iron atoms and for both modifications (see 
Table [U. As the MS simulation shows, the presence of a 
large quadrupole interaction is concealed under the Curie 
temperature by superposition of few subspectra. 

Observing the change in the carbon position upon an- 
nealing the plastically-deformed cementite is hardly pos- 
sible using standard simple approaches to the MS. A first 
sign of existence of the octahedral modification in sam- 
ples and its disappearance with annealing is a change of 
the width of lines in the MS and of the relative distances 
between lines (compare the MS's marked as 3 in Fig. [3]). 

So, to construct an initial model of the local atomic 
structure for a mathematical processing of MS, a simple 
approach does not serve: it is not possible to finally receive 
an adequate fit of the MS. A correct description can be 
obtained using the results of first-principles calculations. 
Certainly, one must bear in mind that first-principles cal- 
culations do not give very accurate magnitude of hyperfine 
parameters and give correctly the tendencies in their be- 
havior and the relations between different sites. 
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